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ABSTRACT 

Aims. This work pursues the investigation of a previously proposed correlation between chemical properties and physical evolutionary 
stage of isolated low-mass star-forming regions. In the past, the A'nh^/A'ccs abundance ratio was suggested to be a potentially useful 
indicator for the evolutionary stage of cloud cores. We aim to study its applicability for isolated Bok globules. 

Methods. A sample of 42 Bok globules with and without signs of current star formation was searched for CCS(2i-lo) emission, the 
observations were complemented with NH3 measurements available in the literature and own observations. The abundance ratio of 
both molecules is discussed with respect to the evolutionary stage of the objects and in the context of chemical models. 
Results. The A'nh^/A'ccs ratio could be assessed for 18 Bok globules and is found to be moderately high and roughly similar across all 
evolutionary stages from starless and prestellar cores towards internally heated cores harbouring protostars of Class 0, Class I or later. 
Bok globules with extremely high CCS abundance analogous to carbon-chain producing regions in dark cloud cores are not found. 
The observed range of A'nhi/A'ccs hints towards a relatively evolved chemical state of all observed Bok globules. 

Key words. ISM: clouds ~ stars: formation - ISM: molecules - radio lines: ISM 



1. Introduction 

Bok globules, named in honor of the astronomer Bart Bok who 
drew attention to those objects and their po ssible role in the 
star formation process (lBok&Reillvl(T947h . appear as small 
and isolated dark clouds. Although the majority of stars in the 
Galaxy is formed in Giant Molecular Cloud complexes, the 
small globules have been recognized as particularly interest- 
ing targets for a study of low-mass star formation, since they 
represent a less complex environment and are therefore more 
easily described by theoretical models. The conditions inside 
Bok globules and the properties of the young stellar objects 
(YSOs) embedded therein have been studied at various wave- 
lengths in several surveys in the past. Therma l radio emis 



sion arising from dust is observed at ce ntimeter dMoreira et al 



1997, [19991) and millimeter wavelengths (Launhardt & HenningI 
19971 hereafter LH97; Henning & Launhardt 1998); submil- 



limeter emission g ives hints for deeply embedded protostars 
(iHuard et al.ll 19991) . With near-infrared o bservations, candidates 
for YSOs are identified and examined (lYun&Clemensl[T99l 
lAlves & Yuiil Il994t iRacca et all |2009|) . Molecular line obser- 
vations provide information about the physical conditions, e.g. 
temperature, density, and magnetic fields, and trace the velocity 
structur e of the clotids. Th ey can be used to identify collapsing 
clouds (IWang et al.lll"995h and reveal the presence of powerful 
molecular outflows dYun & Clemenslll992i 1 199 4b'). A compre- 
hensive overview of star forma tion in different s t ages i s achieved 
by multi-wavelength studies (iLaunhardt et al.l 1201 Oi hereafter 
L2010). 



* Based on observations obtained with the 100-m telescope of the 
MPIfR (Max-Planck-Institut fiir Radioastronomie) at Effelsberg and 
the 64-m Parkes radio telescope. The Parkes radio telescope is part 
of the Australia Telescope National Facility which is funded by the 
Commonwealth of Australia for operation as a National Facility man- 
aged by CSIRO. 



In globules and cloud cores a large variety of molecules 
has been identified. The complex chemical processes leading 
to their formation and destruction are ruled by the local condi- 
tions, including surface reactions on ice-covered grains in the 
coldest and densest parts of the cloud, and are influenced by 
the heating of a forming star within the core or the ionizing in- 
terstellar radi ation close to the cloud surface (review by, e.g., 
Ivan Dishoeck & Blake 1998). While ammonia was the first poly- 
atomic molecu le to be discovered in the interstellar medium 
(review e.g. bv lHo & Townesl[l983|). the CC S radical has been 
identified only two decades ago (Sait o et al.l[T987h and the main 
pathways leading t o its format i on are stiU unde r discussion 
iSuzuki et al.lll992l: lPetriel[T9%l: ISakai et"ani2007l) . although it 
is freque ntly obse rved along with other carbon-chain molecules 
(e.g.. Hirotaetal.l l2009). 

From a study of CCS and NH3 among other molecules, 
ISuzuki et aLld 19921) found that carbon-chain molecules are abun- 
dant in the early evolutionary phases of dark cloud cores, while 
ammonia tends to be more abundant in more evolved, actively 
star-forming regions. This result was also supported by their 
chemical model calculations, where CCS is formed but also de- 
stroyed rapidly, while in turn a replenishment of the molecule 
is impeded by the increasing lock-up of carbon atoms in CO. In 
contrast, NH3 forms in a sequence of very slow reactions and 
reaches its highest abundance in the late evolutionary phases. 
Therefore they proposed the NunJNccs abundance ratio to be 
a possibly useful indicator for the evolutionary stage of star- 
forming clouds. In the following years, the abundance of both 
molecules in a number of Bok glob ules and around Herbig 
Ae/Be stars (IScappini & Codenalll996) . as well as in a larger 
sample of intermediate-mass and low-mass star-forming regions 
with H2O maser emission ( de Gregorio- Monsalvo et aP l2006l 
hereafter dGM06) was examined. However, both studies suffered 
from a small number of CCS detections, which hindered a defini- 
tive judgement about this hypothesis. Focusing on statistics on a 
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Table 1. Overview of the observed globules. 



Object 


RA 


Dec 


Associated 


Distance 




(B1950) 


(B1950) 


IRAS Source 


(Ref.) 




( h m s ) 


(° ' " ) 




[pc] 


CB3 


00 25 59.0 


+56 25 32 


00259+5625 


2500 (2) 


CB6 


00 46 34.4 


-^50 28 25 


00465+5028 


600 


(3) 


CB 12 


01 35 24.3 


+64 47 59 


01354+6447 


800 


(2) 


CB 17 


04 00 32.9 


+56 47 52 


(04005+5647) 


250 


(3) 


CB 22 


04 37 23.5 


+29 49 17 




140 


(2) 


CB 23 


04 40 20.5 


+29 33 26 




140 


(2) 


CB 26 


04 55 56.1 


+52 00 17 


04559+5200 


140 


(3) 


CB 28 


05 03 51.3 


-04 00 18 


05038-0400 


450 


(3) 


CB 34 


05 44 05.7 


+20 59 30 


05440+2059 


1500 


(3) 


CB44 


06 04 43.0 


+ 19 28 19 


(06047+1923) 


700 


(3) 


CB 68 


16 54 27.2 


-16 04 48 


16544-1604 


160 


(3) 


CB 125 


18 12 39.0 


-18 12 14 


(18121-1813) 


200 


(2) 


CB 179 


19 01 57.0 


-05 25 35 


(19018-0525) 


200 


(2) 


CB 188 


19 17 54.0 


+ 11 30 10 


19179+1129 


300 


(3) 


CB 222 


20 32 49.9 


+63 52 00 


20328+6351 


400 


(3) 


CB 224 


20 35 30.9 


+63 42 48 


20355+6343 


400 


(3) 


CB 230 


21 16 50.8 


+68 04 52 


21169+6804 


400 


(3) 


CB 232 


21 35 14.4 


+43 07 17 


21352+4307 


600 


(3) 


CB243 


23 22 51.9 


+63 20 1 1 


23228+6320 


700 


(3) 


CB244 


23 23 48.9 


+74 01 10 


23238+7401 


200 


(3) 


CB246 


23 54 12.0 


+58 17 27 




140 


(3) 


BHR12 


08 07 39.0 


-35 55 54 


08076-3556 


400 


(1) 


BHR13 


08 26 44.1 


-33 36 31 


08267-3336 


400 


(1) 


BHR15 


07 14 27.5 


-43 52 26 


07144-4352 


400 


(1) 


BHR23-1 


08 33 42.6 


-40 28 02 


08337-4028 


500 


(1) 


BHR28 


07 26 20.0 


-50 58 18 




400 


(1) 


BHR34 


08 25 03.4 


-50 30 34 


08250-5030 


200 


(3) 


BHR36 


08 24 15.9 


-50 50 44 


08242-5050 


400 


(1) 


BHR41 


08 26 11.5 


-51 00 39 


08261-5100 


400 


(1) 


BHR55 


09 44 57.0 


-50 52 06 


09449-5052 


300 


(1) 


BHR59 


11 05 03.0 


-61 49 36 




200 


(4) 


BHR71 


11 59 03.1 


-64 52 11 


11590-6452 


200 


(3) 


BHR74 


12 19 21.0 


-66 10 30 




175 


(4) 


BHR86 


13 03 41.4 


-76 44 03 


13036-7644 


180 


(3) 


BHR 1 1 1 


15 38 34.0 


-52 38 21 




250 


(4) 


BHR 121-1 


16 54 59.6 


-50 30 58 


16549-5030 


125 


(4) 


BHR 121-2 


16 55 27.8 


-50 31 01 


16554-5031 


125 


(4) 


BHR 137 


17 18 08.6 


-44 06 17 


(17181-4405) 


700 


(3) 


BHR 138 


17 15 54.0 


-43 24 04 


17159-4324 


225 


(4) 


BHR 139-1 


17 17 15.9 


-43 16 54 


17172-4316 


225 


(4) 


BHR 140-1 


17 19 18.8 


-43 19 24 


17193-4319 


225 


(4) 


BHR 148 


17 01 09.2 


-36 13 59 


17011-3613 


175 


(4) 



Refer ences. (1) Bourke et al. (1995b); (2) iLaunhardt & Hennind 
| |1997|): (3) p- .aunhard t et al., (,20 10,) and private communication; 
(4) lRaccaetairfc009h . 



large number of d ense cores wi t hin the Perseus Molecular cloud, 
a recent work by iFoster et"an (l2009h provided evidence for a 
lower fractional abundance of CCS in protostellar cores com- 
pared to starless cores. 

In this work we are focusing on small Bok globules, typically 
forming only single or few low-mass stars, since they represent 
the most secluded star-forming environments available for inves- 
tigation. An interference by other nearby, possibly high-mass, 
forming stars can be excluded. Moreover, a mixing with ambi- 
ent material which could modify the chemical composition and 
exacerbate a comparison with theoretical models, as it is possi- 
ble in the case of a dense core inside a larger cloud complex, is 
implausible for the isolated globules. An overview of the sam- 
ple examined in this work is given in Sect. |2] followed by a 
description of the observational results and analysis in Sect. H) 



The discussion of the derived physical parameters is presented 
in Sect.|5] concluding remarks are given in Sect.|6l 



2. Description of the sample 

Our sample consists of 42 Bok globules from the catalogs of 
Clemens & Barvainis (1988) for the northern hemisphere and 
Bo urke et al.l (I1995ai) for the southern globules. The northern 
sources have been searched for NH3( 1 , 1 ) and (2,2) emission with 
the Ef felsberg 100-m telescope in the survey of iLemme et al.l 
( 1996, hereafter L96), while the southern globules are included 
in the NH3 survey of Bourke et al. ( 1995b, hereafter BHR95b), 
so that a uniform set of ammonia measurements is available 
for comparison with our data. Table [T] lists the positions of the 
CCS measurement, the IRAS point sources associated with the 
globules (sources in brackets were located outside the telescope 
be am), an d their distance. Following the classification scheme 
of lL20IOl we sorted the objects into th ree groups according to 
the pr esence and evolutionary class (cf. lLadallI987l : lAndre et al.l 
Il993h of embedded YSOs: 

- Group -I: Starless or prestellar globules, or only very young 

and very low-luminosity source embedded. 

- Group 0: Globules harbouring a Class protostar. 

- Group I: Globules containing YSOs of Class I or later. 

The classification of the individual clouds is indicated in Table|2] 
Where too sparse observations of the specific globule prevent a 
reliable identification of the evolutionary group, a presumptive 
classification is given in brackets. It should be emphasized that 
especially prestellar cores or Class protostars will most prob- 
ably remain unidentified when there is a lack of detailed maps 
of molecular emission and/or millimeter continuum emission. A 
detailed discussion for most globules of the sample can be found 
in L2010 and the references therein, while a number of southern 
globules has been classified in iRacca et al.l ( 120091) . In the fol- 
lowing we comment therefore only shortly on the globules not 
described there. 

CB3. A compact submillimeter source without infrared 
coun terpart is located at the center of a bipolar molecular out- 
flow (lYun & Clemensl[r994bl) a nd was sugge sted to consist of an 
aggregate of Class sources bv lHuardetal](l2bOQ) . The source 
IRAS 00259+5625 is located 15" away ^ its near-infrared coun- 
terpart was classified as Class II object (Yun & Clemens 1994^ 
1 1995). T he globule also harbours a H20-maser (Scappini et al 



ll991lidGM 06): it is believed to be rather an evolved interme- 
diate to h igh-mass star-forming region ( ILaunhardt et al.lFl998bt 
ICodella & Bachiller 1999). 

CB 12. An associated IRAS point sourc e is detected at 60 and 
lOOjum, but neither a molecular outflow (.Yun & Clemensl[l992h 
nor dust continuum emission at 1 .3 mm (lLH97h point towards 
the presence of an embedded protostar. 

CB22 and CB 23. Both globules are not detected in the far- 
infrared (no IRAS sources) or in the 1.3 mm dust continuum 
(lLH97h . According to the shape of spectral lines, CB 22 seems 
to be qu iescent while CB 23 sh ows signatures of a possible infall 
motion (" Lee et alJ[T999ll2004l) . 

CB28. Only emission longward of 60jum was de- 
tected by IRAS toward s this globule; searches for outflows 
d^iL&Clemens '1994b') and YSO candidates in the near- 
infrared (Yun & Clemens 1994a) resulted in no detections. 

CB34. Numerous studies suggest multiple star formation 
in the globule. Five submillimeter sources were detected by 
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iHuard et alJ (l2000l) and classified as probable Class ob- 
jects. The aggregate of protostars is associated with a sys- 
tem of multiple outflows (Codella & Scappini 2003) and jets 
dMoreira & YunI Il995h . Evidence for the presence of more 
evolved Class I and Class II objects is provided by the detec- 
tion o f several very red objects in the near-infrared (lAlves & YunI 
119951) . 

CB44. The a s sociate d IRAS source cited by 
IClemens & BarvainisI (Il988l) resides at the very edge of 
the globule. Two 3.6 cm continuum sources without infrared 
counterparts, candidate Class objects, were detected by 
iMoreira et al] 11999) but no detection was possible at shorter 
wavelengths (e.g. at 1.3 mm by LH97). 

CB125. Several IRAS point sources are located in the 
vicinity, but none within the area of highest extinction visi- 
ble on optical images of the globule. Both IRAS 18127-1803 
and IRAS 18122-1818 have rehable fluxes only at 12 and 
25jur n, they were consi dered as candidate pre-main sequence 
star (ILH97I) and YSO (iLee & MversI [T999h . respectively. The 
100//m detection towards IRAS 18130-1824 was considered as 
cirrus em ission due to the lack of submillimeter continuum 
emission dHuard et al.l ll999). Line observations in CS towards 
IRAS 1 8 1 26- 1 820 at the souther n edge of CB 125 sh ow unsuspi- 
cious gaussian line profiles (Lau nhardt et al.llT998ah . Due to the 
the presence of candidate YSOs close to the cloud, we assume 
an evolved stage. 

CB 179. The association with a cold infrared source detected 
only at 60 a nd 100//m, the absence of 1.3 mm continuum emis- 
sion (lLH97h and narrow, weak CO lines (IClemens etalj|1991h 
suggest an early evolutionary stage. 

CB222. Neither is the associated IRAS source detected at 
wavelengths shortward of 60 //m, nor is 1.3 mm dust emission 
found towards this globule (lLH97l) . Together with unsuspicious 
gaussian line profiles measured in CS jLaunhardt et al. 1998a), 
we assume an early evolutionary stage for CB 222. 

BHR 13. The narrow CO line profiles observed towards this 
cometary globule are supposed to arise from cold, quiescent gas 
(jOtrupcek et al."2000), while the associated IRAS source was 
identified as a T Tauri star (Sahu & Sahu 1992). 

BHR 15. The IRAS source associated with this cometary 
globule is detected only at 60 and 100/zm, and CO line 
profiles typical of c old, quiescent gas have been observed 
( IOtrupceket"al]|2000l) . Therefore, an early evolutionary stage is 
assumed. 

BHR 23. ISantos et alj ([1998) conclude from near-infrared 
observations that an aggregate of several YSOs might be har- 
boured b y this globule. High-velocit y wings in CO li nes de- 
tected by lJiwniSretil](|2009) and Otrupcek et al. (2000) might 
indicate the presence of a molecular outflow. Methanol maser 
emission, which is believed to be associated with hi gh-mass star 
formation, was also observed towards this source dWalsh et alj 
[T997I) . 

BHR 28, BHR 59, BHR 74 and BHR 111. These globules do 
not harbour any IRAS sources. Profiles of CO lines observed to- 
wards BHR 74 and the cometary glo bule BHR 28 suggest the 
presence of cold and quiescent gas dOtrupcek et al.l ITOO O*). In 
BHR 59 and BHR 111 the detection of line wings might indicate 
the prese nce of several b lended velocity components or outflow 
motions dOtrupcek et alJl200Q) . However, no more detailed ob- 
servations a re available for thos e globules. 

BHR 41 . ISantos et all d 1 998h classified the two near-infrared 
sources separated by 4" seen at the position of the IRAS source 
as likely Class I objects. Ho wever, no 1.3 mm dust emission was 
detected towards this source dL2010h . 



BHR 137. ISantos et alJ dl998h detected three infrared 
sources at the IRAS position, two of which are likely Class II 
objects. However, the IRAS source is located at the rim of the 
globule and it is unsure if is relate d. Emission in the 1 .3 mm dus t 
continuum has been detected bv lHenning & La unhardll d 19981) . 
as well as a blue wing in the CO emission. IL20 10 suggested a 
rather early evolutionary stage. 

Some of the globules are known to contain multiple sources 
in diffe rent evoluti onary stages, namely CB 3 dYun & ClemensI 
19941 iHuardet al . 2000), CB 34, CB224, CB232, CB243, 
CB244 and BHR 12 (L2010). In such cases, where there are 
sources of different evolutionary stage within the beam of the 
NH3 and CCS observations, we assigned the globule to the evo- 
lutionary group that agrees with the source from which most of 
the mm dust emission arises, i.e., which has the largest reser- 
voir of cold gas and is therefore dominating the detected line 
emission. However, it cannot be excluded that one of the other 
sources also contributes to the detected line emission. Also, the 
maturity of individual globules might be underestimated, be- 
cause the presence of an evolved YSO nearby a pre/protostellar 
core possibly affects its chemistry, especially since in most cases 
the projected separation between the differently evolved sources 
is only a few thousand AU. We have noted these cases in Table|2] 

3. Observations 

Observations in C C S(2i- lo) of all objects from the 
Clemens & B arvainisI dl988l) catalog listed in Table [T| as 
well as NH3(1,1) and (2,2) observations for three globules, 
have been carried out with the Effelsberg 100-m telescope of 
the Max-Planck-Institut fiir Radioastronomie (MPIfR) during 
March 14 and 15, 1999. A maser receiver in the primary focus 
was used with a 1024 channel autocorrelator. The system 
temperature was in the range of 80 and 1 10 K. For the CCS(2i- 
lo) fine at 22.344 GHz the total bandwidth of 6.25 MHz of 
the spectrometer was used, providing a velocity resolution of 
0.08 kms"'. The autocorrelator was split into 2x512 chan- 
nels for the observations of the NH3(1,1) and (2,2) lines at 
23.694 GHz and 23.722 GHz, respectively, corresponding to 
a velocity resolution of 0.15 kms"'. At the used frequencies, 
the beam diameter (FWHP) was 40" . Frequency switching was 
used with a total integration time per source of 30-60 min for 
the CCS line and 5-30 min for the NH3 lines. The pointing 
was checked every 2 hours by cross-scans on nearby quasars. 
Typically, the pointing accuracy was better than 5" at similar 
elevations. The focus was checked on strong continuum sources 
at the beginning of the observing run and after sunset and sun- 
rise. The data were calibrated using cross scans on continuum 
sources with known flux density (Ott et al. 1994). As a primary 
calibrator, NGC7027 was used. The antenna temperatures 
were converted to a main beam brightness temperature Tmb 
by correcting with the elevation dependent gain, the beam and 
aperture efficency. The calibration is estimated to be accurate 
within -15%. 

The 21 sources from the catalog of iBourke et al.l dl995al) 
have been observed with the Parkes 64m telescope during March 
2010 in the CCS(2i-lo) transition, and ten of them additionally 
in NH3(1,1) - namely those sources, where the core or YSO of 
interest was not included in the beam of BHR95b. The 13MM 
receiver and a digital filterbank with a bandwidth of 8 MHz and 
2048 channels were employed, resulting in a velocity resolution 
of 0.05 kms"'. The beam diameter was 1' and the pointing was 
estimated to be accurate within 20". Except of BHR 15, BHR 23 
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Fig. 1. Spectra for the globules with detected CCS(2i- 
lo) emission. For each source a velocity window of 
ca. 20 km s ' around the line is displayed. The derived line 
parameters are given in Table |2] 



-5 5 

velocity [km s ] 



-5 _ 5 

velocity [km s ] 



and BHR28 which were observed in the position switching 
mode, frequency switching with a throw of 1 MHz for the CCS 
and 0.3 MHz for the NH3 line, respectively, was used with inte- 
gration times of 20-60 min. Typical system temperatures ranged 
from 70 to 130 K. After correction of the elevation-dependent 
gain, we applied relative calibration by means of repeated obser- 
vations of BHR71 compared to the measurements of iBHR95bl . 
Hereby, it is ensured that the scale of our measurements and 
those of BHR95b, to be used for calculation of the N^Hj/Nqcs 
ratio, match. Although with this approach the absolute accuracy 
of our calibration is expected to be not optimal, comparison of 
the CCS observations for CB 28, which have been carried out in 
Parkes and EfFelsberg, show consistency. 

The spectra for all globules showing (possible) detections in 
the CCS line are displayed in Fig.[T] 



4. Results 

The results of the observations are summarized in Tables |2] and 
[3] Out of the 42 globules observed, thirteen were detected in 
the CCS(2i-lo) line with a signal-to-noise ratio of at least three. 



Additional five globules show emission at a slightly lower signif- 
icance level, but we believe them to be real detections since their 
LSR velocities are in agreement with those of ammonia lines ob- 
served by BHR95b. In contrast, NH3(1,1) remained undetected 
only towards three globules of our sample. 

From the further discussion we exclude CB 3, CB 34 and 
BHR 23 because they likely contain clusters or high-mass star- 
forming regions, as well as BHR 13, BHR 59 and BHR 1 1 1 be- 
cause of their very uncertain evolutionary stage, and BHR 12 be- 
cause of the presence of two equally massive sources in different 
evolutionary groups (the source is however included in Figs. |2] 
and|6]l. Considering only the globules with fairly reliable evolu- 
tionary group (i.e. group given without brackets in Table|2|l, CCS 
emission is present in 70% (seven out of ten objects) of the glob- 
ules belonging to the starless or prestellar globules of group -I, 
in 40% of the objects of group (clouds containing Class pro- 
tostars; four detections) and 33% of the group I objects (YSOs of 
Class I or later embedded; three out of nine detected). The over- 
all detection rate (43%) of CCS(2i-lo) in our sample is there- 
fore larger than found in some earlier studies of low-mass star- 
forming regions, globules and dark clouds (about 18-33% in the 
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Table 2. Derived CCS(2i-lo) line parameters and column densities for CCS and ammonia, and resulting abundance ratio NnhJNqcs- 



Source 


Group 


t^LSR 


Tmb 




.1 


du 




A'ccs 


A'nh, (Ref.) 


A'nhj/A'ccs 






[km s 


J 


[mK] 


[kms-'] 


[mK km s ' ] 


[kms-'] 


[10'-cm-2] 


[lO'^cm-2] 




CB3 









<96 














< 1.42 


1.54 ±0.26* (1) 


>104 


CB6 


I 






<90 














< 1.32 


0.29 ± 0.07* 


>22 


CB 12 


-I 






< 117 














<1.70 


1.45 ±0.24* (1) 


>86 


CB 17 


-I" 


-4.65 d 


-0.02 


233 ± 37 


4.S -1- 
u.tj in 




1 12 + 


12 


0.44 + 


u.uo 


4.89 ±0.53 


13.84 ±3.74 


201 ±47 


CB 22 


-I 


-1.03 d 


-0.01 


233 ± 37 






36 -1 


- 7 


U.UO It 


u.uo 


1.42 ±0.29 


3.64 ±0.87(1) 


257 ± 80 


CB 23 


-I 


4.86 ± 


0.01 


375 ± 39 


0.17 + 


0.02 


66 ^ 


- 8 


0.1 1 + 


0.03 


2.39 ±0.28 


10.1 ± 1.1 (1) 


422 ± 67 


CB 26 


I 






< 114 














< 1.67 


0.48±0.18* (1) 


>28 


CB 28 


-I 


9.01 ± 


0.02 


223 ± 40 


0.29 ± 


0.04 


69 il 


-9 


0.26 ± 


0.05 


3.48 ±0.45 


1.95± 1.18(1) 


56 ±34 


CB 34 





0.17 ± 


0.05 


119 ±37 


0.60 ± 


0.14 


57 ± 


10 


0.59 ± 


0.15 


2.64 ±0.53 


0.92 ±0.44* (1) 


34 ±17 


CB44 


(-1) 






< 132 














<1.95 


0.65±0.15* (1) 


>34 


CB 68 









< 120 














< 1.48 


12.7 ±0.9(1) 


>860 


CB 125 


(I) 






< 132 














< 1.88 


6.01 ±1.63(1) 


>320 


CB 179 


-I 






< 102 














< 1.49 


0.59±0.15* (1) 


>40 


CB188 


I 


6.14± 


0.04 


111±35 


0.44 ± 


0.08 


52 ± 


10 


0.42 ± 


0.09 


2.21 ±0.41 


1.10±0.33* (1) 


50 ±17 


CB 222 


-I 


-2.70 d 


:0.05 


102 ± 25 


0.34 ± 


0.12 


36 ± 


12 


0.31 ± 


0.14 


1.54 ±0.53 


0.65 ±0.08* 


42± 15 


CB 224 


0* 






< 144 














<2.61 


6.23± 1.17(1) 


>239 


CB 230 




2.41 ± 


0.03 


190 ±49 


0.33 ± 


0.12 


66 ± 


16 


0.30 ± 


0.13 


2.79 ±0.66 


0.71 ±0.24* (1) 


26±9 


CB 232 


-I'' 


12.50 ± 


0.03 


138±41 


0.23 + 


0.07 


34 ^ 


- 9 


0.19 + 


0.09 


1.42 ±0.38 


0.54±0.21* (1) 


38 ±15 


CB243 








< 114 














< 1.67 


2.77 ±1.04(1) 


>166 


CB244 





4.18 ± 


0.03 


162±41 


0.31 ± 


0.07 


54 ± 


10 


0.29 ± 


0.07 


1.96 ±0.37 


8.08 ±2.37** (1) 


412 ± 143 


CB246 


-I 


-0.65 d 


lO.Ol 


462 ± 39 


0.24 ± 


0.03 


120 ± 


11 


0.21 ± 


0.03 


5.67 ±0.49 


7.70± 1.01 (1) 


136±18 


BHR12 


I-hO^ 


6.36 ± 


0.07 


190 ±68 


0.35 ± 


0.15 


71 ± 


23 


0.33 ± 


0.16 


2.75 ±0.88 


5.4 (2) 


196 ±62 


BHR13 


(I) 






<294 














<7.14 


0.68 ± 0.22* 


>10 


BHR15 


(-1) 






<221 














<5.40 


0.7 (2) 


>13 


BHR23 


(I) 






<206 














<5.01 


4.07 ± 0.67* (2) 


>81 


BHR28 


(-1) 






<209 














<5.08 


3.2 (2) 


>63 


BHR34 









<244 














<5.95 


0.91 ±0.50* 


>15 


BHR36 


I 


5.23 ± 


0.11 


125 ±61 


0.94 ± 


0.17 


124 ± 


26 


0.93 ± 


0.17 


4.13 ±0.86 


7.22 ±2.01 


175 ± 60 


BHR41 


IS 






<198 














<5.07 


1.02 ±0.54* 


>20 


BHR55 


0" 


-5.36 d 


- 0.05 


173 ±93 


U.Z. / It 


U. iU 


iUU It 




U.Z.J It 


1 1 

U. i 1 


3.38±2.11 


5.30 ±1.45 


157 ± 107 


BHR59 


(0) 






<221 














<5.37 


2. 19 ±0.32* (2) 


>41 


BHR71 


0' 


-4.52 d 


-0.08 


209 ± 78 


U.o / ± 


n 1 zL 


1 -t- 




U.oO It 


1 zL 


6.63 ±1.15 


12.25 ±2.75 


185±51 


BHR74 


(-1) 


J .07 II 


- U.U7 


1 1 1 ± 45 


0.52 ± 


0.22 


61 ± 


18 


0.50 ± 


0.22 


2.19 ±0.66 


1.3 (2) 


59±18 


BHR86 









<264 














<7.23 


5.0 (2) 


>69 


BHRlll 


(0) 


-0.07 d 


:0.03 


186 ±44 


0.30 ± 


0.06 


59 ± 


12 


0.28 ± 


0.07 


2.42 ± 0.48 


3.6 (2) 


149 ± 29 


BHR 121-1 


I 






<212 














<5.14 


0.42 ±0.23* 


>8 


BHR 121-2 


I 






<201 














<4.90 


0.48 ±0.35* 


>10 


BHR 137 


(0) 






<245 














<5.95 


4.3 (2) 


>72 


BHR 138 









<223 














<5.41 


< 0.46* 




BHR 139-1 


I 






<211 














<5.12 


<0.45* 




BHR 140-1 





-0.63 d 


:0.07 


221 ±55 


1.21 ± 


0.16 


285 ± 


31 


1.21 ± 


0.16 


11.1 ±1.20 


20.4 (2) 


184 ±20 


BHR 148 









<282 














<3.52 


<0.80* 





Notes. (*) Optical thin approximation used to derive A'nh3, (**) A'nh, derived from available map and map-averaged T^ot and Tex, since NH3 peak 
lies outside the CCS beam. Main core is prestellar (SMMl+2), but there is a low-luminosity Class I source 20" away, slightly outside the beam. 

Main core is Class (SMMl), but there is a Class I source 30" away (not inside the telescope beam). Contains two Class I IR sources 10" 
apart. Main core is most likely prestellar (SMM), but there is a Class I source 12" away. '•'^^ Main core is most likely prestellar (SMMl), but 
there is a Class I source 10" away. A Class and Class I YSO with cores of similar mass separated by 20" inside the beam. Two close (4" 
separation) infrared sources of Class I dSantos et al. 1998). '''' Two close lines with comparable width a;0.27 km s"' in CCS. Double protostellar 
core with two IR sources 17" apart. 

References. (D iLemme etaP dl996h : (2) lBourke et"an jl995br) . 



studie s of IScappini & CodellallT996l ldGM06l and ISuzuki etal] 
II992I) . and closer to the 50% detection rate of a survey of dense 
cores in the Perse us Molecular Cloud dRosolowskv et al.ll2008 l: 
[Foster et"al]|2009h . For ammonia, strong emission in the sense 
that the hyperfine structure of the NH3(1,1) transition is detected, 
is present in 80% of the objects from group -I, and in 70% and 
22% of the objects of group and I, respectively. 



4.1. Analysis of CCS lines 

In Table |2] Column (1) lists the object number from the 
Clemens & Barvainis (1988) or Bourke et al. (1995a) catalog, 
(2) the assigned evolutionary group (cf. Sect. |2]i, (3) the LSR- 
velocity, (4) the main beam brightness temperature, (5) the ob- 
served line width (FWHM) including instrumental broadening 
(of 0.08 kms ' and 0.05 kms ' for the CB and BHR sources, 
respectively), and in (6) the integrated intensity of the line. 
Column (7) contains the nonthermal linewidths, (8) the calcu- 
lated total column densities for CCS, (9) the ammonia column 
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Table 3. Results of the NH3(1,1) and NH3(2,2) observations. 



Source 


t'LSRCl.l) 




At'obsd.l) 


!^LSR(2 


.2) 




AL'obs(2,2) 




Tex 




[km s-'] 


[K] 


[kms-'] 


[km s 




[K] 


[kms-1] 


[K] 


[K] 


CB6 


-12.45 ±0.09 


0.15 ±0.05 


0.62 ±0.1 8 






<0.14 






6* 


CB 17 


-4.65 ±0.01 


1.71 ±0.10 


0.34 ±0.02 


-4.70 ± 


0.06 


0.17 ±0.04 


0.79 ±0.1 8 


2.8 ±0.6 


5.4 ±0.4 


CB 222 


-2.71 ±0.04 


0.29 ±0.05 


0.72 ±0.10 






<0.11 






6* 


BHR13 


6.44 ± 0.06 


0.30 ±0.07 


1.25±0.18 












6* 


BHR34 


4.78 ±0.04 


0.50 ±0.10 


0.69 ±0.11 












6* 


BHR36 


5.27 ±0.02 


1.21 ±0.07 


0.52 ±0.04 










1.9±0.5 


4.3 ±0.2 


BHR41 


5.46 ±0.05 


0.36 ±0.08 


0.92 ±0.16 












6* 


BHR55 


-5.50 ±0.01 


1.45±0.11 


0.51 ±0.04 










1.2±0.3 


5.1 ±0.4 


BHR71 


-4.39 ±0.02 


3.56±0.11 


0.71 ±0.03 










1.5±0.3 


8.0 ±0.6 


BHR 121-1 


-12.61 ±0.04 


0.26 ±0.06 


0.61 ±0.08 












6* 


BHR 121-2 


-12.52 ±0.05 


0.27 ±0.07 


0.68 ±0.15 












6* 


BHR 138 




<0.21 














6* 


BHR 139-1 




<0.20 














6* 


BHR 148 




<0.36 














6* 



Notes. (*) Tex is adopted as average value from lLemme et al.] ( 119961) and lBourke et alj ( Il995bh . 



densities from own observations or collected from the litera- 
ture, and (10) the abundance ratio A^nHj/A^ccs- Source velocity 
i/LSR, hnewidth Ai/obs, peak intensity T^b and integrated inten- 
sity r rmbdf were derived by Gaussian fits of the lines using 
CLASSQ For the error of r„b, the root mean square (rms) noise 
of the spectrum, and for the other quantities the standard devia- 
tions of the Gaussian fits to the lines are listed in Table |2l The 
average rms noise level of the obtained spectra is 38 mK for 
the CB sources and 72 mK for the BHR sources. As upper lim- 
its for the non-detections, a peak intensity of 3 rms is given. In 
the case of CB 34 and BHR 55 the observed CCS line was not 
well approximated by a Gaussian profile, therefore the value in 
Table |2] represents the intensity integrated under the actual line 
rather than under the fitted profile. 

The total column densities A^ccs were calculated under the 
assumption of a local thermodynamic equilibrium (LTE): 



gj 



Mw:rzY)i ^^^^ 



(1) 



where v, j denotes the frequency of the transition, go and gj are 
the statistical weights of the ground state and of the upper ro- 
tational level, respectively, Ej is the energy of the upper level, 
Q(Ts^) denotes the partition function for an excitation tempera- 
ture Tex, Aji is the Einstein coefficient for spontaneous emission 
and Tv the optical depth. We assumed the emission to be opti- 
cally thin, thus the relation between optical depth and main beam 
brightness temperature r^b from the radiative transfer equation 
is given by 



E^Jk = 1.61 K (IWolkovitchetal.lll997h . and Q^A5 K)=23.80 
(iLai & Crutcheiil2000h . the column density ensues as follows: 



A^ccs [cm ^] 



2.92 X 10' 



r^bdw [K km s 



(3) 



A 4 K higher excitation temperature would result in a mod- 
erate increase of about 40% for the column density. For the 
not detected lines the upper limit for J TmbdL' was estimated as 

1 .06rn,bAL', where Av was taken to be the average of the detected 
line widths, 0.32 kms ' for the CB sources and 0.64 kms ' for 
the BHR sources. 

4.2. Ammonia lines 

For most of the sources in our sample, ammonia column den- 
sities were derived in the papers of L96 and BHR95b. For our 
own ammonia observations, column densities were derived as 
outlined in the following. If hyperfine structure for the NH3(1,1) 
transition was sufficiently strong, a fit with the available proce- 
dure within the CLASS package was performed to obtain opti- 
cal depth of the main line Tm(i,i) and intrinsic Hnewidth Ai/im. 
Otherwise a Gaussian profile was fitted to the main line. The de- 
rived line parameters are summarized in Table[3] where Auobsci.i) 
refers to the intrinsic hnewidth for CB 17, BHR 36, BHR 55 and 
BHR 71, and to the width of a single Gaussian fitted to the main 
component otherwise. Further analysis was p erformed follow- 
ing standard proced ures and assumptions (e.g. lUngerechts et aH 
119801 : IStutzkilll984h . The excitation temperature Tgx was calcu- 
lated from 



1 



.hVijlkT„ 



1 



(2) 



where Tbg - 2.73 K is the cosmic background temperature. 
Since only one rotational transition of CCS was observed, the 
excitation temperature could not be determined. We adopted 
Tex ~5 K, equal to the average value fo r a number of dark cloud 
cores as found by ISuzuki et alj (Il992h . Restricting to the mea- 
sured CCS(2i-lo) transition with A2.-1 



' part of the program package GILDAS (Grenoble Image and Line 
Data Analysis Software), see http://www.iram.f r/IRAMFR/GILDAS/, . 



mb(l,l) 



/b (1 - 



The 



(4) 



with the beam filling factor f\, and brightness temperature of the 
main line rmb(i,i)- The column density in the (LI) levels then 
follows as 



(1,1) - 



(5) 



with the rate coefficient An = L67 x 10"^ s ' and frequency 
vii of the NH3(1,1) transition; denotes the line strength at 
the main component, the value of which depends on the intrinsic 
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linewidth of the hyperfine components blended in the main line. 
The complete ammonia column density is then estimated by 



NH3 



1 



IIK 5 



14 



(6) 



Since we did not observe the NH3(2,2) transition in most 
cases, we could deriv e the rotational temperature (see e.g. 
lUngerechts et"anil98(]l) only for CB 17 (r^t = 8.1 ± 0.5 K). For 
the other sources we adopted a rotational temperature of 10 K 
for normal and 13 K for cometary globules (BHR 12, BH R 13, 
BHR15, BHR 28), i.e. average values from the papers oflL96l 
and BHR95b, for the calculation. 

Where the hyperfine structure was weak or only the main 
line of the NH3(1,1) transition was detectable, so that the optical 
depth could not be derived, the column densities were estimated 
in the optically thin approximation with 



16;rv3 
^(U) = — T— 



1 



I + ghvnIkTc 



hvu fhAu 1 



mb(l,l) 



i)dv, (7) 



where J 7'inb(i,i)dL' is the intensity integrated over the main com- 
ponent of the NH3(1,1) transition. For the excitation tem perature 
in equation |7] the average value found by L96 and lBHR95bl 
Tex = 6 K has been assumed. The total column density can then 
be calculated from equation |6] adopting Tiot - 10 or 13 K as 
mentioned above. Upper limits were estimated with the same ap- 
proach as described for the CCS lines, with an average linewidth 
of the NH3(1,1) main component of 0.85 kms"'. This optical 
thin approximation has also been applied for sources which were 
detecte d, bu t for which no column density was derived in the pa- 
pers oflL96l and'BHR95b' Note that for CB 230 and CB 232 the 
NH3 measurement positions of lL96l are o ffset by about one beam 
from the submillimeter cores detected b v ' L20 1 Ol (identical to the 
CCS measurement position), and hence A^nh, can be expected to 
underestimate the column density towards the cores. 

For the calculation of the column densities given in Table |2] 
the beam filling factors /b were adopted in the following way: 
for sources included in L96 and BHR95b the source sizes de- 
rived there were used to calculate the beam filling factor to be 
applied to our Parkes CCS measurements, assuming that CCS 
and ammonia emission have comparable spatial extents (which 
is an approximation, see e.g. Codella & Scappini 1998a). For the 
ot her so urc es, we use d the average source sizes from the papers 
oflL96land lBHR95bl to calculate an average beam filling factor 
of /b = 0.81 for the BHR, and 0.69 for the CB objects. 

5. Discussion 

5. 1 . Column densities and fractional abundances 

The column densities of CCS and NH3, shown in Table |2] do 
not exhibit strong dependence on the evolutionary group of the 
corresponding source, although the highest values for A^ccs are 
found for group objects. Column densities of NH3 versus those 
of CCS are displayed in Fig. |2]( without the omitted sources men- 
tioned in Sect. H)). The lower right part of the plot is occupied 
exclusively by sources undetected in CCS (i.e. only upper limits 
for A^ccs); all sources with A^nh, < 1.5 x lO'^^cm"^ have am- 
monia column densities obtained with the optical thin approx- 
imation. Sources with reliable detection of both molecules are 
found preferentially in the upper left side of the plot, as if the 
highest CCS column densities are present (or excited) only in 
globules with high ammonia column density, while lower val- 
ues of A^ccs can be found in objects with a wider range of A^nh, ■ 



However, no clear distinction is visible between the distribution 
of the presumably younger objects of group -I and the actively 
star-forming globules (groups and I). 

The large scatter of the column densities Ukely results from 
differences in the hydrogen densities of the individual clouds, 
therefore a comparison of fractional abundances is required for a 
reasonable discussion of possible evolutionary trends. Hydrogen 
column densities or masses for several globules were derived 
from dust emission bv lL20ld[ from C'^O measur ements bv lL96l 
'Wang et al.' (1995) and 'Vilas-Boas et alj (11994 . and from ex- 
tinction by Racca et al. (2009). Some globules of our sample are 
included simultaneousl y in thre e of these studies; for these the 
column densities A^h of lL201Cil (derived by dividing the hydro- 
gen masses from their work by the proton mass and the phys- 
ical area enclosed by the 1.3 mm dust emission in the maps) 
and L96 are found to be systematically larger than those of 
IWang et al] (Il995h by a factor of 7.6 and 2 (within 20%), re- 
spectively. Thus, we assume the same conversion factors for the 
rest of th e globules i n the ac cording works to bring them on the 
scale of Wang et al. (1995), but we do not use th is approach 
for Racca et al. (2009) and Vilas-Boas et all (11994 since there 
is an overlap of only one source with the other studies. For 15 
globules with reliable evolutionary stage this yields fractional 
abundances Nccs/Nh of 7 x 10"" to 6 x 10""' and NnhJNh be- 
tween 2 X 10"'' and 1 x 10"^. For CCS, the fractional abundance 
averaged over globules of the same evolutionary group is very 
similar for group -I, and I (about 3 x 10""^). For NH3 it is 
highest for group (6 x 10"^) and lowest for group I (1 x 10"^), 
but this variation should not be considered significant due to the 
large scatter between individual globules and their small num- 
ber. These estimated fractional abundances are comparable with 
those of several chemical models from the literature and an own 
model, which will be described later in Sect. 15.21 at evolutionary 
times close to 10^ yrs. Comparing with the fractional abundances 
observed for dense cores in cloud complexes from Foster et al.l 
(2009), A^NH,/A^H from the globules is similar, but Nccs/Nn is 
about one order of magnitude lower than in the dense cores. 
However, as pointed out earlier in this section, the hydrogen col- 
umn densities of the globules vary up to one order of magnitude 
between different studies and therefore the fractional abundances 
of CCS and NH3 estimated here have to be considered uncertain 
to the same degree. 



5.2. NI-I3/CCS ratio in comparison with chemical models 

In contrast, the relative abundance A'nhs/A^ccs can be discussed 
independent of the hydrogen densities under the assumption that 
both molecules trace similar cloud regions, which is likely be- 
cause of comparable excitation conditions and a similar tele- 
scope beam for both observations. Figure [3] shows A^nh,/A^ccs 
versus evolutionary group. Altogether, we derive abundance ra- 
tios from about 20 to 860, while the observable range is limited 
to ca. A^NH,/A^ccs < 2000 by our CCS detection hmit. 

Compared to samples of earlier papers, the A^nhj/A^ccs ratios 
derived for isolated Bok globules are on average similar to those 
of dense cores in the Perseus molecular cloud (Rosolowskv et al] 
2008), but a factor of two hi gher than those in the dark clouds 
studied bv lSuzuki et al.l (1199 2). This corresponds to the fact that 
despite a comparab le range of ammonia column densities in the 
ISuzuki et alJ (Il992h sample and our globules, we do not find ex- 
tremely high CCS c olumn densities typ ical of carbon-chain pro- 
ducing regions (e.g. iHirota et alj2"009l where carbon-chain pro- 
ducing regions are defined as having A^nhj/A^ccs ^ 10). 
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N(CCS) [cm ] 

Fig. 2. Column densities of NH3 versus CCS. Arrows indicate 
upper limits; filled and open circles represent objects with evi- 
dence for ongoing star formation (group and I) and no signs of 
current star formation (group -I), respectively. Grey symbols de- 
note sources with uncertain evolutionary group (given in brack- 
ets in Table |2]). 



The number of sources is too small to allow for a detailed sta- 
tistical analysis, but taking into account only sources with reli- 
able evolutionary stage, the following picture arises: within each 
evolutionary group, a relatively large range of abundance ratios 
A^NH,/M:cs spanning about one order of magnitude is observed, 
but there is comparatively little variation of the values between 
the different groups. Fig.[3]possibly suggests a slightly decreas- 
ing tendency of A^nhj/A^ccs going from group to group I glob- 
ules, but this trend is only marginal. Thus, it can be concluded 
that the ratio Nj^hJNccs is rather similar across the Bok glob- 
ules observed in this study, despite them harbouring YSOs in 
different evolutionary stages. 

Specifically, the A^NH,/Atcs ratio does not increase going 
from the presumably youngest objects of group -I to the most 
evolved sources of group I. This finding does not easily fit into 
the general picture of most chemical models, according to which 
CCS as an early-phase molecule is expected to decrease in abun- 
dance rapidly around 10*' yr, while the slowly forming ammonia 
reaches its maximum abundance in a later stage of the chemi- 
cal evolution - leading to the anticipation of steadily increasing 
NuhJNccs ratio along with the evolutionary stage of the glob- 
ules. For comparison, we show in Fig. |4] the Nnh,/Nccs ratio 
with respect to the chemical age of a cloud, calculated from the 
evolution of CCS and NH3 abundances from four chemical mod- 
els in the literature: Suzuki et alj (119921 dash-dotted line) and 
IScappini et al.l (11998 . dashed line) use pseudo-time-dependent 



chemical models in which the g as density is c onstant within 
time. In contrast, iBerginI ( l2000l solid line) and ' Aikawa et alj 
(1200 iL dotted line) take into account the dynamics of a collaps- 
ing core, as well as depletion of species from the gas phase 
onto, and desorption from, dust grains. However, both models 
do not include reactions on grain surf aces except of H2 forma- 
tion and ion-electron recombination in I Aikawa et al.l (12001). As 
initial conditions all models take hydrogen in molecular form 
and carbon and sulphur as ions, the initial abundances are taken 
as those typical of diffu se cloud s or depleted by a certain factor 
from solar abundances; IBerginI 12000) allow the cloud to evolve 
for 1 .5 X 10^ yr at constant density and take the chemical compo- 
sition after this time as initial values for the dynamically evolv- 
ing core. All models consider the chemical evolution in a region 
at a constant temperature of 10 K and shielded from external UV 



CO 

o 
o 



X 

z 




gr. —I gr. gr. I 

evolutionary group 

Fig. 3. Abundance ratio A^nh,/A^ccs versus evolutionary group; 
arrows indicate lower limits and grey symbols sources with un- 
certain evolutionary group (slight horizontal offsets around the 
group positions are solely for better visibility of individual data- 
points). 



10° r 



10' 



B 10" r 



z 



10 



10' r 




10' 



lO'^ 



10* 
time [yr] 



10° 



10° 



Fig. 4. NH3/CCS ra tio from ch emical 
(see Sect. 15.2b of Suzu ki et al 



model calculations 
11992, dash-dotted Hne), 



ScaiDiDini et al. ( 1998, dashed line),|Bergin (2000, solid line) and 
lAikawa et al.l (120011 dotted line), and this work (10 K warm 
group -I model marked with open circles, 15 K group model 
with filled circles and 25 K group I model with asterisks. The 
observed range of A^nHi/A^ccs is designated by the shaded area. 



radiation. lAikawa et al.l (1200 1*) follow the chemical evolution of 
an infalling fluid element in a collapse according to the Larson- 
Penston solution (we show he re their result f or a collapse slowed 
down by a factor of 10). From lBerginI (|2000|) we show the model 
for a collapse with ambipolar diffusion and dust grains covered 
with a CO mantle. Both papers examine several variations of 
their models (e.g. grain properties, collapse timescales), which 
result in somewhat different time evolutions of the molecular 
abundances, but do not deviate significantly from the examples 
shown in Fig.|4] 

Despite the different approaches, all models agree upon a 
rapid increase of the A^nh,/A^ccs ratio by three orders of magni- 
tude at an evolutionary time of several 10^ yr, caused by a fast 
decrease of the CCS abundance due to depletion onto grains, 
destruction by reactions and missi ng replenishme n t. The first 
broad peak seen in the results of IScappini et al.l (Il998h and 
lAikawa et alJ (l200ll) is due to a very slow formation and in- 
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Fig. 5. Upper part: Evolution of the NH3 fractional abundance 
resulting from our model described in Sect. 15.21 for group -I 
model (10 K, solid line), group model (15 K, dashed line) and 
group I model (25 K, dotted line). Lower part: The same for 
CCS. 



crease of the NH3 abundance, while CCS is still being formed 
efficiently. 

In general, a comparison of the absolute values observed 
with those from modelling has to be considered with caution. 
The initial conditions assumed and the exact starting point for 
reactions defined in chemical models might not necessarily be 
in good agreement with real globules, and certain scatter has 
to be expected due to the natural fluctuations of initial con- 
ditions in the variety of globules. Nevertheless, comparing the 
measured values for A^nhs/A^cs abundance ratio with the mod- 
els suggests that the observed objects might all be in an evolved 
state matching the region after 10^ yr despite their different evo- 
lutionary stages. In addition, we employed an own chemical 

model to calculate abundances of NH3 and CCS for sources 
of the three groups -I, and I. For this purpose a recent 
chemical network inc luding a variety of grain surface reactions 
(ISemenov et al ]|2010^ was used together with a simplified phys- 
ical model. Representative physical parameters for each group 
were taken from the globule study bv lL2010l i.e. the density is 
described by a power law function of radius (with exponent be- 
tween -1 and -2) in the outer region and flattening towards the 
core center. We used a temperature of 10 K for group -I, 15 K 
for group 0, and 25 K for group I sources, respectively, based 
on observations and models of the mass-averaged dust tempera- 
ture (iL2010i and references therein). Starting from "low metal'" 



atomic abundances (iLee etalJll998h as initial values, the evolu- 
tion of NH3 and CCS abundance relative to H was calculated 
separately for the typical physical conditions of each group. We 
refrained from constructing a piecewise model with sharp tran- 
sitions from one group to the next at defined timesteps, since 
in reality a smooth transition of physical conditions can be ex- 
pected, and the lifetimes of the different groups are not tightly 
constrained. Instead, we discuss qualitative differences of the 
three models. 

The resulting evolution of fractional NH3 and CCS abun- 
dances is displayed in Fig.|5] In our chemical model, CCS forms 
fast in a neutral-neutral reaction of CCH and S, and reaches 
its maximum concentration early. NH3 is formed partly on dust 
grains through surface hydrogenation of nitrogen, and partly in 
the gas phase mainly via a sequence of ion-molecule reactions. 
Desorption from grain surfaces occurs mainly through cosmic- 
ray heating and cosmic -ray-induced far-UV photons in the cold 
and shielded environments considered here. After a few 10^ yrs 
CCS is efficiently removed from the gas phase by freeze out 
onto dust grains or destruction by reactions with oxygen, and 
rapidly drops in abundance. For the 15 and 25 K warm group 
and group I models, prior to this late-time depletion of CCS a 
number of competitive surface reactions become active, which 
enhance the CCS (re)production and lower the NH3 abundance, 
resulting in a temporary drop of the NH3 abundance at about 
10^ yrs and a slightly increased CCS peak abundance compared 
to the 10 K model of group -I sources. For other carbon-chain 
molecules, a similar enhancement in moderately warm environ- 
ments ('warm carbon-chain chemistry') has been pointed out re- 
cently (e.g. Sakai et allllOOSl: lAikawa et al.' 20081: iHassel et al.l 
[2008). These differences in chemical evolution of our models 
originate primarily from the different temperatures, while the 
difference in density profiles has smaller influence. For simplic- 
ity, we will refer to the models of the three groups (-1, and I) 
by their temperature (10, 15 and 25 K) in the following. 

As a result, in the early phase the ratio A^nhj/A^ccs evolves 
very slowly from values close to unity to few dozen for all three 
models. For the 10 K model, the abundance ratio starts to in- 
crease faster around 10^ yrs and crosses the range observed in 
the Bok globules of our sample within 2.5 x 10^ yrs. For the 
somewhat warmer 15 K model, A^nh,/A^ccs increases earlier and 
more rapidly, while for the 25 K model the abundance ratio re- 
mains low (< 10) until > 1 X 10^ yrs before increasing rapidly 
to > 10^ like the other two models. In Fig. |4] the evolution of 
Nj^hJNccs for the three parameter sets of our chemical model 
is designated by open and filled circles (10 K and 15 K, respec- 
tively) and asterisks (25 K). 

Comparing with the observational data shown in Fig. [3] the 
observed range of Ni^hJNccs ratios could be compatible with 
the rapidly increasing parts of our three model curves at evolu- 
tionary times of few 10'*-10'' yrs. However, the 10 K and 15 K 
model are close to the lower Nt-mJNccs limit derived from our 
observations for most of their evolutionary time, while the 25 K 
model remains at significantly lower values up to 10^ yrs. Thus, 
if one assumes that the actual chemical ages of Bok globules 
may possess a certain scatter, one could also expect to find glob- 
ules with very low A^nHi/A^ccs ratios (< 5) if the group I model is 
applicable - but such low ratios are not confirmed in our sample. 
The range of fractional abundances estimated from the observa- 
tions in Sect. 15. II agrees for NH3 with both 10 and 15 K model 
for times > 10"^ yrs, while the NH3 abundance in the 25 K model 
remains slightly below the observed range for the whole time 
span covered by the model. For CCS, the peak abundances of 10 
and 15 K model around 10'*-10^ yrs fall into the observed range. 
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while the 25 K model exceeds it. Although having to keep in 
mind the uncertainties of the estimated fractional abundances, 
this may be another indication that the 25 K model is not in 
agreement with the observations, while 10 and 15 K warm model 
fit them equally well. 

Originally, the choice of 10-25 K for th e models of the dif- 
ferent groups wa s motivated by mode ls (e.g.. lShirlev et alj2002l: 



iGalh et al 



ps wa s motivatea by moae is (e.g.. laniriev et auzuu^t 
2002 : Bergin et"aLll2006h and observations (e.g. for 



CB 244 in Stutz et al] 2010h of the dust temperature. In general, 
dust temperatures are expected to be elevated (~15 K) close to 
the cloud surface heated by the interstellar radiation field, and 
lower (~10 K) in the moderately dense envelope where CO line 
radiation cools the gas efficiently. In the center of prestellar cores 
they may be as low as 5 K (depending on the external radi- 
ation field), or accordingly warmer in the case of an internal 
heating source (20-30 K in IShirlev et al]|2002l) . In the dense 
(hh > 10^ cm"^) inner regions gas and dust are expected to be 
well coupled through collisions and hence similar in tempera- 
ture. 

In contrast, the beam- (BHR95b) or map-averaged (lL96h 
NH3 rotational temperatures (which represent a good estimate of 
the kinetic te mperatur e at the low temperatures prevalent in glob- 
ules, see e.g. Stutzki 1984) are 9-16 K for globules of all evo- 
lutionary stages (L96.; .BHR95b.). This most likely results from 
the fact that the warm cores of Class and more evolved ob- 
jects have typical sizes of few thousand AU (L2010'), and hence 
comprise only a small fraction of the NH3 beam area for the av- 
erage globule distance of our sample. In addition, the innermost 
warmest and densest parts of the core may not be traced well by 
NH3 due to high optical depth. Thus, with the large beam of the 
NH3 observations considered here, cool (~ 10-15 K) gas from 
the moderately dense envelope ca n be expected to dominate the 
signal (see also lCrapsi et al.ll2007l where the central gas temper- 
ature drop in a prestellar core is detectable with interferometric, 
but not single-dish observations). 

Moreover, a scatter of temperatures in the range 10-15 K for 
the bulk of gas may also result from differences in size, density 
and local UV background radiation for the i ndividual globules , 
and not only from their evolutionary stage. In'F oster et al.! (2009) 
the kinetic temperatures from NH3 of protostellar and starless 
cores cover a similar range as those for our globules, and the 
slight differences between the star-forming and starless group 
detected by them may be discernable only for sample sizes sig- 
nificantly larger than ours. 

If indeed, as suggested by our chemical model, small tem- 
perature differences of ~ 5 K could already perceptibly af- 
fect the progression of the A^nh,/A^ccs ratio around 10^ yrs, 
then the superposition of A^nhJA^ccs evolution curves for sev- 
eral initial temperatures could lead to a relatively large scatter 
of A^NH,/M:cs ratios within one evolutionary group, while si- 
multaneously smearing out differences between the evolutionary 
groups (differences of other initial physical parameters as den- 
sity may also contribute, but we do not address them here in 
detail since the observational data is too sparse for comparison). 
This would agree with the observation that the (beam-averaged) 
A^NH,/A^ccs ratio is rather similar across globules of all evolu- 
tionary stages and exhibits a relatively large scatter within each 
group. Additional contributions to the scatter of the Nt-iaJNccs 
abundance ratio among globules of a single group may arise 
from the significant fraction of globules containing adjacent ob- 
jects of different evolutionary stages (Sect. |2]l, indicative of se- 
quential star formation, and a range of ages even among the ob- 
jects compiled in one evolutionary group. 
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Fig. 6. Abundance ratio A^nhs/A^ccs versus bolometric tempera- 
ture of embedded sources. Arrows indicate lower limits of the 
abundance ratio, grey symbols mark Tboi of combined spectral 
energy distributions (cf. Sect. 15.3b . 



To test this hypothesis in detail, it would be necessary to dis- 
entangle the influence of gas temperature and age (i.e. evolu- 
tionary group), for which our dataset is significantly too small. 
It could also be of interest to evaluate a possible increase of gas 
temperature and change of chemistry on small scales close to 
forming protostars, for which observations at much higher reso- 
lution are needed. 

5.3. Abundance ratio and bolometric temperature 

In the previous years it has been pointed out that the evolution 
of young stellar objects is well traced by the bolometric temper- 
ature Tboi, defined as the temperature of a blackbody with the 
same mean freque n cy as the spect r um of t he observed source 
(iLadd et al.1 119911: iMvers & Laddl Il993h . IChen et al.1 dl995h 
demonstrated a tight correlation between Tboi and the age of pro- 
tostars and pre-main-sequence stars in particular for embedded 
sources. For a number of young stellar objects in globules of our 
sa mple, the bo lometric temperature has been evaluated bv L2010l 
and lRacca et a l. (2009). However, Tboi is mainly useful to track 
the evolution of Class and later sources, since for most of the 
coldest, presumably youngest sources (e.g. CB246, BHR 137) 
Tboi is not known due to a spectral energy distribution not well 
constrained at shorter wavelengths. 

We show our measured A^nh VA^ccs ab undance ratio versu s 
Tboi adopted from the works by L2010' and iRacca et al ] (I2009h . 
in Fig. |6] In the case of CB 232 and CB 243 (as well as BHR 12 
denoted in grey colour), where both a prestellar core and a 
Class I source are included in the beam, Tboi represents the value 
for the combined spectral energy distribution of the adjacent 
sources, which is in both cases dominated by the Class I source 
while the molecular emission is expected to emerge mainly from 
the prestellar cores. However, in both globules the separation of 
prestellar core and Class I source equals the typical size of a 
prestellar core (8000 AU in L2010); in this case the presence of 
an evolved YSO might influence the chemistry of the prestellar 
core and accordingly both may be not completely unrelated. For 
a definite judgement about this question, a mapping of the fine 
emission at small scales would be required. 

From most chemical models in the literature as shown in 
Fig- III ™ increase of the A^nh,/A^ccs ratio towards warmer bolo- 
metric temperatures, which in turn are thought to represent in- 
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creasingly older and more evolved sources, would be expected. 
Despite the scarceness of our dataset, there are no indications of 
such an increase. Instead, the NnhJNqcs abundance ratio shows 
a tentatively trend of decrease towards increasing bolometric 
temperature. However, due to the small number of globules with 
both known bolometric temperature and detected CCS emission, 
an interpretation is even more difficult than in Sect. 15.21 

5.4. Nonthermal linewidths 

The observed linewidths A^obs include an instrumental broad- 
ening of Afinstr = 0.08 kms"' (Effelsberg) and 0.05 kms"' 
(Parkes) which can be subtracted under the assumption of a 
Gaussian profile for the line as well as for the spectrometer: 

Av'^Avl^-Avl, . (8) 

In the same way, the thermal contribution to the linewidth, 

/8hi2l7\~ 

AWtherm = \ , (9) 

V m 

with m denoting the molecular mass of the molecule, was sub- 
tracted to obtain the nonthermal linewidths Ai /nth in Ta bled We 
used the rotational temperature Tiot derived bv lL96l and lBHR95bl 
from their ammonia observations as an estimate for the kinetic 
temperature Tyn. Where such values for the individual globules 
were not available, an average temperature of Tkin = 10 K was 
assumed. The typical thermal contribution for the CCS and NH3 
linewidths is 0.09 and 0.16 kms respectively. It should be 
pointed out that, because the main component of the NH3(1,1) 
line consists of several usually blended hyperfine components 
distributed over 0.7 kms linewidths derived from a single 
Gaussian fit to the main component are dominated by the sep- 
aration of the blended hyperfine components rather than the in- 
trinsic linewidth for small linewidths. An intrinsic linewidth can 
be in principle estimated from such fits, but suffers from large 
uncertainties for the typical small linewidths observed in Bok 
globules. Therefore we discuss nonthermal linewidths of only 
those globules with sufficiently well detected main and satellite 
lines of NH3(1,1), for which intrinsic linewidths can be derived 
directly from fits of the hyperfine structure. 

For both NH3 and CCS, the majority of nonthermal 
linewidths is found in the range 0.2-0.5 kms ' for all evolu- 
tionary groups. In addition, large nonthermal linewidths 0.6- 
1.2 kms ' are observed in globules of group and I, while very 
small linewidths <0.1 kms ' are present exclusively in glob- 
ules of group -I. The nonthermal linewidth can be understood 
as a measure for turbulence of the medium and the existence of 
unresolved velocity components within the area covered by the 
telescope beam. An increased linewidth Ai/nth for actively star- 
forming regions could be expected due to turbulence and out- 
flow motions induced by the presence of a protostar. However, 
for globules at larger distances a physical region of proportion- 
ally larger dimension is enclosed by the beam, and could there- 
fore result in larger linewidths. We are fairly confident that our 
results are not biased by this effect, since no clear relationship 
between linewidth and distance is discernable in our sample. 

Among the starless or prestellar sources of group -I, the 
globules with extremely small nonthermal linewidths (CB22 
and CB 23) are distinguished from the rest by an absence of 
IRAS and millimeter sources, i.e. there is no evidence for any 
embedded sources. In contrast, three of the group -I globules 
with intermediate Av^h harbour more evolved Class I sources 



in the vicinity of the prestellar core, while the remaining glob- 
ules are poorly studied. The group and I sources with large 
nonthermal linewidths exhibit not very Gaussian-shaped profiles 
(BHR 140-1, BHR36) or hints of at least two velocity compo- 
nents (BHR71, also BHR 55). However, they do not differ dis- 
tinctly from the rest of the group and I globules in properties 
like presence of outflows, multiplicity, etc. 

In the majority of globules, the nonthermal CCS linewidths 
are smaller than those of NH3, but also the opposite case or 
similar linewidths are observed. A possible explanation could 
be that both molecules trace not e xactly similar cloud regions , 
as in the cases demonstrated by ICodella & Scappiiiil ( Il998bl) 
and Lai & Crutcher (200 0), where the most dense parts of cores 
were preferentially detected in NH3, while CCS emission arises 
from a more extended surrounding region. Such central de- 
pletion holes in evolved sources have also been observed for 
other carbon-bearin g species, e.g. for th e CCH radical (e.g., 
.Beuther et aL ,20081 : iPadovani et al. I l2009h . which resembles in 
many properties, and is thought to be a precursor of, CCS. Since 
both NH3 and CCS possess a roughly comparable critical den- 
sity (~ lO^'-lO"* cm"-*), a difference in the emitting region is most 
likely a result of different spatial distributions. 

This in turn might mean that both molecules trace re- 
gions with somewhat different physical conditions, and also re- 
sult in an incorrect estimate of the thermal component for the 
linewidths when assuming the same kinetic temperature for both 
CCS and NH3. Altogether, it is difficult to draw firm conclusions 
regarding the origin of the spread of nonthermal linewidths due 
to the sparse dataset and availability of only single-point mea- 
surements. 

6. Summary 

Observations of 42 Bok globules in CCS and partially in NH3, 
supplemented with NH3 measurements from the literature (lL96l 
BHR95b) enabled us to assess the abundance ratio of both 
molecules for 18 objects (and derive lower limits for additional 
21 globules) in different evolutionary stages ranging from star- 
less and prestellar globules (designated as group -I) over clouds 
containing Class YSOs (group 0) to globules harbouring young 
stellar objects of Class I and later (group I). In the following we 
summarize our main results: 

1. The CCS(2i-lo) line is detected in 18 of 42 Bok globules. 
The detection rate is highest in the objects of group -I (70%) 
and decreases towards the globules of group and I (40 and 
33%, respectively). Ammonia is detected (together with lit- 
erature data) in 39 globules. 

2. We derive Nj^hJNccs ratios in the range 26-422, plus one 
lower limit of 860, for isolated Bok globules. In particular, 
within the Umits of our survey, we find neither extremely 
low abundance ratios (< 10) like those typical of carbon- 
chain produ cing regions in dark cloud (Suzuki et al. 199^ 
iHirota etldj r2()09l) . nor extremely high abundance ratios 
(several 10^) as expected for evolved sources from chemi- 
cal models for cloud cores by different authors (Fig.|4]i. 

3. We do not observe an increase of the A^nhj/A^ccs abundance 
ratio going from the starless and prestellar globules towards 
evolved globules containing Class I or later sources, nor from 
lower towards higher bolometric temperatures of the em- 
bedded YSOs, as would be expected from various chemical 
models from the literature and this work (Sect. lS!2b . Instead, 
the ratio exhibits a considerable scatter but is roughly con- 
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Slant across all evolutionary groups (Fig. O, with a ten- 
tatively - however, statistically not significant - decreasing 
trend from globules containing Class protostars (group 0) 
towards globules with Class I or later YSOs (group I). 

4. An own chemical model (Sect. I5.2l i indicates that a slight 
temperature increase (15 instead of 10 K) could affect the 
NuhJNccs ratio perceptibly at evolutionary times around 
10^ yrs. Since NH3 rotational temperatures vary in the same 
range, this may suggest that the observed roughly constant 
but scattered distribution of (beam-averaged) A^nHi/A^ccs ra- 
tios could result from a superposition of evolutionary tracks 
for different initial globule temperatures. In contrast, a 25 K 
warm model seems less likely, since neither such high gas 
temperatures, NH3 and CCS fractional abundances in agree- 
ment with the model, nor very low Nt-mJNccs ratios (< 5), 
as expected from the model even at late evolutionary times, 
are observed. 

5. The smallest nonthermal linewidths of CCS (~ 0.1 kms '), 
indicating a very low level of turbulence, are detected only 
in two globules of group -1 without any associated infrared 
or millimeter sources. In contrast, broad CCS lines with 
evidence for multiple velocity components are found only 
among star-forming globules (groups and 1). Beyond that, 
no firm conclusions about the nonthermal linewidths can be 
drawn due to the small sample size. 

We conclude that our observed NnhJNqcs abundance ra- 
tios derived from single-dish observations with relatively large 
beam, although related to the evolutionary state of the embed- 
ded objects, cannot be alone and straightforward interpreted as 
an evolutionary tracer for isolated Bok globules. Another major 
problem hampering this study is the immediate vicinity of ob- 
jects in different evolutionary stages encountered in many of the 
globules (see also L2010). which makes it difficult to assess the 
amount of emission arising from each source or a possible mu- 
tual influence on physical and chemical conditions. In addition, 
the assumption of NH3 and CCS tracing the same spatial regions 
is an approximation and may not hold for all globules. The distri- 
bution of CCS - central depletion hole or enhanced in the warm 
core region of evolved sources - is of particular interest. Well- 
resolved mapping of the line emission of both molecules for a 
large sample of globules, assessment of the gas temperature, and 
a position-dependent evaluation of the Nt^nJNccs ratio, will be 
required to address the mentioned aspects in full detail. 
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